Zero finite-temperature Drude weight within the ID half-filled Hubbard model 
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The Hubbard model is a problem of wide physical interest that possibly is the most studied lattice 
model of correlated electrons. In spite that in one dimension (ID) it is solvable by the Bethe ansatz, 
at half filling its finite-temperature T > transport properties remain poorly understood. Here we 
combine that solution with symmetry to show that within that prominent T = ID insulator the 
charge Drude weight D(T) vanishes for T > in the thermodynamic limit. This result is exact and 
clarifies a long-standing open problem. It rules out that at half filling it is an ideal conductor in the 
thermodynamic limit. Evidence is provided that it is a normal resistor. 

PACS numbers: 71.10.Fd, 72.80.Sk, 72.20.-i, 71.10.Hf 
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The nature of the exotic transport properties of one- 
dimensional (ID) correlated systems at finite tempera- 
ture has been a problem of long-standing interest [UtM- 
The ID Hubbard model is solvable using the Bethe ansatz 
(BA) [5-7]. This technique has been useful in the calcu- 
lation of static properties. However, it has been difficult 
to apply to the study of transport at finite temperature. 

The real part of the charge conductivity as a function 
of the energy fiw and temperature T has the general form, 



fictitious flux 4>/N a , the model Hamiltonian reads, 



a{u,T) = 2ttD(T)S(lu) 



(1) 



Here the charge stiffness D(T) characterizes the response 
to a static field, within linear response theory. Moreover, 
a reg (uj,T) describes the absorption of light of frequency 
uj. At T > the system can behave as an ideal conductor 
with D(T) > 0, a normal resistor with D(T) — and 
(To = lim^^o OregJ (uj, T) > 0, and an ideal insulator with 
D(T)=a o = Um- 

The solvable ID Hubbard model has D(T) > in the 




This is 
10 1, since 



metallic phase at finite temperature T > 
consistent with a general result due to Mazur 
some of the model's conserved quantities have nonzero 
overlap with the charge current operator. On the other 
hand, the charge transport in its T = insulating phase, 
which corresponds to N = N a electrons and lattice sites 
and thus n = N/N a = 1 electronic density, is not well 
understood. For instance, whether in the thermodynamic 
limit the charge stiffness D(T) vanishes or is finite for 
T > and n = 1 remains an open issue 0, Q • 

Here we fully clarify that unsolved problem by showing 
that in the thermodynamic limit, D(T) = for T > 
at n = 1 and U/t > 0. Here t is the nearest-neighbor 
transfer integral and U the on-site repulsion. Our result 
definitively establishes that the half-filled ID Hubbard 
model is not an ideal conductor. In the presence of a 



H 



h.c. 



N a 



i=i 



(2) 

Here cj a creates an electron of spin projection a at site 

j and hj, a = cj j(T Cj t(T . 

Symmetry plays a major role in our study. Recently it 
was found that for U/t ^ the global symmetry of the 
Hubbard model on a bipartite lattice and thus in ID is 
[SO(4)®U{\)]/Z 2 [nj]. In addition to the 77-spin and spin 
517(2) symmetries in 50(4) = [5*7(2) ® SU(2)]/Z 2 , it 
contains a hidden U(l) symmetry whose generator eigen- 
value can be chosen to be the number of rotated-electron 
doubly plus unoccupied sites [ll|, [l2| , which we denote 
by 25^. The 77-spin (and spin) and 77-spin projection 
(and spin projection) of the ID Hubbard model energy 
eigenstates are denoted by 5^ and 5^ (and S s and S"), 
respectively. The S a and 5^ values of the lowest-weight 
states (LWSs) and highest-weight states (HWSs) of the 
77-spin and spin algebras are such that S a = —5^ and 
S a = 5£, respectively. Here a = 77 for 77-spin and a = s 
for spin. Moreover, 25^ = 2[5 r) + M'] where M' is the 
number of 77-spin-neutral pairs of rotated-electron dou- 
bly and unoccupied sites 12j, which in ID equals the BA 
number M' of Ref. @|. 

The model's BA solution refers to a subspace spanned 
either by the LWSs or HWSs. Here we use the LWS 
representation of that solution and call Bethe states the 
energy eigenstates that are LWSs of both the spin and 
77-spin 5*7(2) algebras. For them the numbers, 
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= S v --(N a -N) = 0,l,...,2S n .. 



- 5, 



-(N f -N i )=0,l,...,2S s 



(3) 



vanish. We denote the energy eigenstates by \^>u & ). A 
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FIG. 1: (a) The degenerate energy spectrum of the S v = 
1; M' — and S n = 0; M' = 1 states and the current spectra 
of (b) the metallic S v = 1; S* = -1; M' = states, (c) half- 
filling S v = 1;S* = 0;M' = and S 1 ,, = = 0; M' = 
1 states, and (d) metallic S v = l',SZ = 1; M' = states 
considered in the text for U/t = 8 and JV a = 90. 



non-LWS energy eigenstate |'4 r j,{ A ) is generated from the 
corresponding Bethe state p ) as follows, 
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I.I". 



(4) 



FIG. 2: The degenerate energy spectrum of the S v = 2; M' = 
states, S v = 1; M' = 1 states, and S v = 0; M' = 2 states 
considered in the text for U/t — 8 and iV a = 30. 



Here, 



= [n a \] JJ[25 a + l-j / " 



a = rj, s . 



(5) 



are normalization constants, 



N„ 



3 J J • 

S"4 c i,t ' 



5+ = 



£4 

3=1 



4. C J:t ' 



(6) 



and /§ a = (S^)^ are the 77-spin (a = 77) and spin (a — 
s) off-diagonal generators, and Ia and l A stand for the 
set of numbers [2S*, S^, S s , n v , n s ] and [2S*, S n , S s , 0, 0], 
respectively. Within our notation, l A refers to limiting 
values of the general index I a such that n n = n s = 0, 
which arc those of the Bethe states. Moreover, I stands 
for all remaining quantum numbers beyond I a needed to 
uniquely define an energy eigenstate |^z,i A ). 

Except for a constant pre-factor, the charge current 
operator J p equals the 77-spin current operator J a * , 



J p 
J 



(e)J; J'i 
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(J j = l 



(1/2) J , 



C j+l,cr C j,CT 



(7) 



where e denotes the electronic charge. The ID Hubbard 
model in the presence of a fictitious flux 4>/N a , Eq. 
is solvable by the BA, which provides the <j>/N a depen- 
dence of the Bethe states, E u o^(cf>/N a ). The Bethe-state 

expectation values of the current operator J of Eq. $7$ , 



I.I" 



dE lil0 J<j>/N a ) 
d(<P/N a ) 



4=0 • 



(8) 



can then be extracted from the BA solution [13|, |14 |. 
From the use of Eq. (HJ) one finds that the expectation 
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values of such a current operator of the non-LWSs of a 
multiplet-7/-spin SU(2) tower of states can be written as, 



-4-3-2-101234 



1 



Jim = 7T ^i^\{s-rJ{s+T^ ul o ) 



(9) 



where the normalization constant C n is that given in Eq. 
([5]) for a = r] and |^,;<^) is the tower ?7-spin LWS. 

The charge stiffness D(T) can be calculated as a ther- 
modynamic quantity [l], Q , 



D(T) 



1 d 2 F 

2N n dtf 



\<t>=o + 



2TN a f^ 



Pl,l & (jl,ll 



(10) 



The first term involves the free energy F and vanishes 
in the present thermodynamic limit. In the second term 
Pi,i& — e~^ El,,& jZ is the usual Boltzmann weight, Z 
stands for the partition function, and Jz,; A is the current 
expectation value given in Eqs. ([5]) and ©. 

The n = 1 states are of two types: (a) S v = S2 = 
Bethe states l^z^) inside the BA solution subspace 
and (b) energy eigenstates outside that subspace with 
S v = 1,2,3,... integer and S* — 0. For simplicity we 
consider the n = 1 energy eigenstates of zero spin den- 
sity yet our results hold as well for the half-filling states of 
finite magnetization. It trivially follows from the fy-spin 
SU(2) symmetry algebra that the 77-spin current expec- 

tation value J, ?„ = ;o I J CT - 1$; ; o ) of the zero-r;-spin 

'A 'A 

S n = SZ = Bcthe states vanishes. The simple rela- 
tion between the 77-spin and charge current operators re- 
ported in Eq. (J7|) then plays a major role in our study. 
It implies that the corresponding charge current and cur- 
rent Ji t i° A of Eq. © vanish for such states. Hence they 

give no contribution to the charge stiffness D(T), Eq. 
(|TU| . Thus whether in the thermodynamic limit D(T) 
is finite or vanishes for T > is fully determined by 
the contributions from the above energy eigenstates with 
S v = 1, 2, 3, ... integer and SZ — 0. 

Such n — 1 states belong to the same 77-spin tower as 
metallic Bethe states l^^) with S v = — S* — 1, 2, 3, .... 
On using the methods of Ref. [3] , one finds by means of 
the BA solution that the current J^i° & , Eq. ©, carried 
by such states is in general finite. Fortunately, one can 
use a suitable operator algebra to exactly evaluate the 
current Jz i ; A , Eq. ©, carried by the corresponding half- 

0. 



filling states with S v 



1, 2, 3, ... and thus S* 



Specifically, the systematically use of the commutators, 



where, 



N a 
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,t 



(11) 



(12) 



and J = ( J + )t allows to relate the current expecta- 
tion value Ji,i A of Eq. © to the expectation value J^o , 
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FIG. 3: The current spectra of (a) the metallic S v = 2; S% = 
-2;M' = states, (b) metallic S n = 2; = -I'M' = 
states, (c) half-filling S v = 2; S% = 0; M' = states, S v = 
= 0;M' = 1 states, and S v = S* = 0; M' = 2 states, 
(d) metallic S 1 ,, = 2; — 1; A/' = states, and (e) metallic 
S n = 2; — 2; M' = states considered in the text for 



U/t 



and iV a = 30. 
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Eq. (|SJ), of the corresponding 77-spin-tower LWS. After a 
suitable operator algebra involving commutator manipu- 
lations one finds, 

Jl,l ix =C{2S ri ,n n )J l ^. (13) 

The coefficient C(l, I) is such that, 

C(lJ) = -C(l,l-l), !<l/2 
C {1,1/2) = for 1/2 integer, (14) 

where that 0(1,1/2) = follows from the first equality 
for 1 = 1/2 and we have denoted 2Sr, and n v by / and /. 

The half-filling energy eigenstates have numbers S ri = 
n v = 1,2,3,..., so that C(2S n ,n v ) = C(2S V ,S V ) = 0, as 
given in Eq. (|14[) . Importantly, it then follows from Eq. 



(fl~3| that in spite of such n = 1 states being generated 
from S v = —S* = 1,2,3,... metallic states, Jy A = 
for all of them. It was found above that also J, ,o = 

'A 

for the remaining half- filling states with S v = S* = 0. 
One then concludes that the charge current expectation 
value vanishes for all n = 1 states. On the other hand, it 
follows from Eq. (TTU1) that D(T) vanishes in the thermo- 
dynamic limit provided that the charge current expecta- 
tion value of all half-filling energy eigenstates vanishes as 
well. Hence we have just showed that at half filling D(T) 
vanishes for T > in the thermodynamic limit. 

Only the coefficient C(2S ri ,n v ) = C(2S V ,S V ) = is 
needed for our study. The general expression of the coef- 
ficient C(2S V , rirj) for n v = 0,1,2,3 and any ?7-spin value 
Sr/ > n v /2 (for n n > 3 it becomes too cumbersome) is: 



C(lJ) 



IILiIjJ - 2* ] - (I - 1)(2[(J + 1) - (I - l)])'- 1 - (1 - SjMl - 2')(1 2)(2[(J + 1) - (I - 2)}) 



1-2 



[ii]nL[f+i-j] 



l = 1,2,3, l > l; l = 2S„, l = n r 



(15) 








1 


2 


3 


4 


5 


6 


7 


1/2 


1 


-1 














1 


1 





-1 












3/2 


1 


1/3 


-1/3 


-1 










2 


1 


1/2 





-1/2 


-1 








5/2 


1 


3/5 


1/5 


-1/5 


-3/5 


-1 






3 


1 


2/3 


1/3 





-1/3 


-2/3 


-1 




7/2 


1 


5/7 


9/16 


1/7 


-1/7 


-9/6 


-5/7 


-1 



TABLE I: The coefficient C(2S V , n v ) on the right-hand side of 
Eq. (|13[) for the ?7-spin-tower states of 77-spin up to S v = 7/2. 
At half filling one has that S v = 0,1,2,3 is an integer and 



0, so that n„ 



S v > 0. Moreover, J t , 



*^7> C (^Srj , Srj) 

,0=0 for 5„ = 



0. 



and J;,; z 
0. 



for 



Combining the expression of Eq. (fT5")) with the relation 
C(l,t) = -C(l,l- I) provided in Eq. (O for I < 1/2 
we have calculated the coefficient C(2S v ,n v ) of all states 
with ?7-spin S v < 7/2 whose values are given in Table |U 
Finally, we have combined our results with the BA 
techniques of Ref. iTT| to derive the currents, Eqs. © 
and ©, of related half- filling and metallic energy eigen- 
states with S n = 0, 1, 2 and 2S* c h = 2,4 for U/t = 8. We 
considered all states with 2S^ = 2[S n + M'] = 2 and 
thus two holes in the c momentum band [151 . That in- 



cludes the three types of S* = 0,±1; S n 



= 1; W = 



states and the S„ 



■,, - -S* = 0; M' = 1 states. The 
energy and current spectra of such states are plotted in 
Fig. [TJ Furthermore, we have considered states with 
2S' r l = 2\S n + M'} = 4 c band holes. That includes the 



five types of S* = 0,±1,±2; S n 

Sr] — 

0; M' = 



= 2; M' = states, 
1; M' = 1 states, and 
2 states with two and 

'1 r,„A t'2 



three types of S* — 0, ±1 
two types of S v = —S* 

one occupied BA quantum number J^ 1 and J£ of Ref. 
0, respectively. The energy spectrum of all such states 
is degenerate and is plotted in Fig. [5J The current spec- 
tra of the energy eigenstates with 2S^ — 2[S V + M'] = 4 
and 77-spin S v = and S v = 2 are plotted in Fig. El The 
spectra of Figs. 1(a), (b) and 2, 3(a) were calculated from 
the BA for N a = 90 and N a = 30, respectively. 

Recently the finite-energy behavior of correlation func- 
tions of ID correlated systems 16-2(| has been found to 
differ signicantly from the linear Luttinger liquid theory 
predictions [2lj. Here we have considered the related 
problem of the exotic T > charge transport proper- 
ties of the half-filled ID Hubbard model. We have shown 
that its charge Drude weight D(T) vanishes for T > in 
the thermodynamic limit. Combining that exact result 
with the finite cto found numerically in Ref. 0] provides 
evidence that in the thermodynamic limit that model is 
for T > and at n — 1 a normal resistor. 
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